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Abstract 

We use chiral perturbation theory to study the extrapolations necessary to make physical pre- 
dictions from lattice QCD data for the electromagnetic form factors of pseudoscalar mesons. We 
focus on the quark mass, momentum, lattice spacing, and volume dependence and apply our re- 
sults to simulations employing mixed actions of Ginsparg- Wilson valence quarks and staggered sea 
quarks. To determine charge radii at quark masses on the lattices currently used, we find that all 
extrapolations except the one to infinite volume make significant contributions to the systematic 
error. 
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I. INTRODUCTION 



As computing resources and numerical algorithms improve, first principles determination 
of hadronic observables will be possible with lattice QCD. In the foreseeable future, these 
numerical determinations will rely on effective field theories to address systematic errors in 
lattice data. There has been considerable effort to understand and compute effects from 
the finite volume of the lattice, the discretization chosen for fermions, and most notably the 
quark mass dependence of observables. To address systematic error related to the treatment 
of the fermionic determinant, quenched chiral perturbation theory 0, 0, HI, and partially 
quenched chiral perturbation theory jE IE IE 0| have been developed. In the respective 
approximation made in the former case, the determinant is replaced by a constant leading to 
uncontrolled systematic error, while in the latter the determinant is computed but with larger 
quark masses used than in the propagators connected to external legs. Only through this 
latter approximation can connection to real QCD observables be made; moreover, effective 
field theory is required to make this connection systematically. 

In this work we study the extraction of the electromagnetic charge radii of pseudoscalar 
mesons from lattice QCD within the framework of partially quenched chiral perturbation 
theory. We consider the various extrapolations in quark mass, momentum, volume and lat- 
tice spacing needed to extract the radii. Experimentally, the pion charge radius is rather 
well determined from pion scattering off atomic electrons [E EE]- As such it can be used 
as a crucial test of lattice and effective field theory methods. Away from small momentum 
transfer, the pion form factor has been experimentally probed from the virtual pion cloud 
of the nucleon however, extraction is limited by model dependence assumed in extrap- 
olating the experimental data to zero virtuality. Ultimately lattice methods will enable first 
principles QCD calculation of meson form factors over a wide range of momentum transfer. 

The original lattice QCD calculations of the pion form factor were pursued by two 
groups [l2|, HH . Since these pioneering calculations there have been various further computa- 
tions and refinements using improved actions, larger volumes, and different lattice fermions; 
for recent investigations see, e.g., [3, EE EE E3|- These calculations are all limited by 
the quenched approximation. There have been, however, recent lattice calculations that 
include dynamical quarks [lEE^- The results with dynamical quarks in [l8j] use a mixed 
lattice action of domain wall valence quarks on staggered sea quarks, and the lightest pion 
mass ~ 300 MeV is arguably within the chiral regime. Indeed such mixed action simulations 
are currently popular due to both the publicly available MILC configurations [2(J, and the 
desirable chiral symmetry properties of domain wall fermions [21] (or more generally, of 
Ginsparg- Wilson fermions 1221 ). For other observables calculated with such mixed actions, 
see B A HE EE H H Ufi HQ. The low-energy effective theory for mixed lattice 
actions is a partially quenched chiral perturbation theory even when the valence and sea 
quark masses are degenerate. An additional reason for us to study the pion form factor 
is to learn about the low-energy constants of mixed action chiral perturbation theory |33| . 
At the order we work, only one new parameter C m i x [see Eq. (J21|) ] enters in the continuum 
extrapolation of the pion charge radius. We find, however, that the data in E3 allow for 
only a rough estimate of the parameter C m i x , a problem that can be remedied with more 
lattice data for differing valence quark masses, or at different lattice spacings. 

The organization of the paper is as follows. First in Sec. |Hj we review the basics or 
partially quenched chiral perturbation theory. In Sec. IIII[ we derive the meson form factors 
at one-loop order in partially quenched chiral perturbation theory. Our result improves upon 
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an earlier calculation [34| by using a computationally judicious choice for the quark charges. 
Next in Sec. IIVI we include the effects from the lattice discretization for the case of a mixed 
action of Ginsparg- Wilson valence quarks and staggered sea quarks. The finite volume 
corrections are presented in Sec. Extrapolations are considered in Sec. I VH where we use 
our results to investigate the chiral, momentum, volume, and continuum extrapolations of 
meson form factors. Lastly we conclude with a brief summary of our work, Sec. IV111 



II. PARTIALLY QUENCHED CHIRAL LAGRANGIAN 

In partially quenched QCD, the quark part of the continuum Lagrangian is written as 

9 

c = J2QAP+^q)Qj- (i) 

3=1 

The nine quarks appear in the vector 

Q = (u, d, s, j, l,r,u, d, 5) T , (2) 



that transforms in the fundamental representation of the graded group SU(6\3) |35|, |3 
The quark mass matrix is given by 

rn Q = diag(m„, m d , m s , mj, m h m r , m u , m d , m s ), (3) 

to maintain the cancellation of path integral determinants from the valence and ghost sectors. 
Effects of dynamical quarks are present due to the contribution of the finite-mass sea quarks. 
Additionally we choose to work in the isospin limit in the valence and sea sectors: m d = m u 
and mi = rrij. 

The light quark electric charge matrix Q is not uniquely defined in partially quenched 
QCD |37(. The only constraint one must impose is that the charge matrix Q has vanishing 
supertrace. Following 3^ 3^|, we use 

Q = diag (q u , q d , q s , q j} q h q r , q u , q d , q s ) , (4) 

along with the condition qj + qi + q r = 0. QCD is recovered in the limit of degenerate valence 
and sea quarks only for the particular choice: q u — qj — |, and q d = q s = qi = q r = — §■ 
Letting the charges be arbitrary, however, enables us to track the flow of charge in loop 
diagrams. 

For massless quarks, the Lagrangian in Eq. (0) exhibits the graded symmetry SU(6\3)l(8) 
SU(6\3)r <S> U{l)v that we assume is spontaneously broken to SU(6\3)v <8> U{l)y. The 
low-energy effective theory of partially quenched QCD is written in terms of the pseudo- 
Goldstone mesons emerging from spontaneous chiral symmetry breaking. At lowest order 
in the chiral expansion, the dynamics of these mesons can be described by the 0(p 2 ) La- 
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C = ^-str (D^D^) - - str (m^S + m Q ^ + ad^ d^ + (5) 



Here p ~ where p is an external momentum. 
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where 

S = exp ) , (6) 

*=( M i), (7) 



X M 

/ = 132 MeV, and we have defined the electromagnetic gauge- covariant derivative -D M S = 
<9 M £ + ieA^Q, £]. The str() denotes a supertrace over flavor indices. The M, M, and x are 
matrices of pseudo-Goldstone bosons with quantum numbers of qq pairs, pseudo-Goldstone 
bosons with quantum numbers of qq pairs, and pseudo-Goldstone fermions with quantum 
numbers of qq pairs, respectively. Upon expanding the Lagrangian in (JSJ) one finds that 
quark basis mesons with quark content QQ' have masses 



rn 



X 

< ><r ~ J2 



— (m Q + m Q/ ). (8) 



The flavor singlet field appearing above is given by $ = str($)/ v / 6. Just as in chiral 
perturbation theory, but in contrast to the quenched case, the singlet field is rendered heavy 
by the strong U(1)a anomaly and is integrated out of the theory. The resulting flavor 
neutral propagators, however, deviate from simple pole forms jsj. We do not display these 
propagators here as they are not explicitly needed in our final results for meson form factors. 

Additionally there are three terms in the 0(p A ) Lagrangian 

C = a 4 str(D M SD At S t )str(mJ ? S + m Q S t ) + a 5 str(£) At SD At S t (mJ ? S + mQS t )) 

-ta 9 str {L^D^Djtf + R^D^DJZ) , (9) 

that contribute to meson form factors at tree level. Here L^, R^ v are the field-strength 
tensors of the external sources, which for an electromagnetic source are given by 

V = R pv = eQ(d^A u - dvAp) + ie 2 Q 2 [A fJi ,A u ]. (10) 

Unlike quenched chiral perturbation theory, where the low-energy constants are distinct 
from those in chiral perturbation theory, the partially quenched parameters in Eq. Q are 
the dimensionless Gasser-Leutwyler parameters of chiral perturbation theory [4(| , which can 
be demonstrated by matching. 



III. FORM FACTORS IN INFINITE VOLUME 

The electromagnetic form factor Gx of an octet meson (f>x is required by Lorentz invari- 
ance and gauge invariance to have the form 

(Mp')\MMp)) = eG x {q 2 )(p + p')^ (11) 

where p (p') is the momentum of the incoming (outgoing) meson, and g M = (p' — p)^ is the 
momentum transfer. Conservation of electric charge protects it from renormalization, hence 
at zero momentum transfer e Gx(0) = Qx, where Qx is the charge of (f> x - The charge radius 
r x is related to the slope of Gx(q 2 ) at q 2 = 0, namely 
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FIG. 1: Loop diagrams contributing to the octet meson charge radii in partially quenched chiral 
perturbation theory. Octet mesons are denoted by a dashed line, singlets (hairpins) by a crossed 
dashed line, and the photon by a wiggly line. 
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FIG. 2: Wavefunction renormalization diagrams in partially quenched chiral perturbation theory. 

Charge conjugation implies the form factor relations: G n +(q 2 ) = —G^-i^q 2 ), Gk+(q 2 ) = 
—G K -(q 2 ), and G K o(q 2 ) = —G^o(q 2 ), as well as G n o(q 2 ) = G v (q 2 ) = 0. 

To calculate the charge radii to lowest order in the chiral expansion one has to include 
operators of £ in Eq. © to one-loop order [see Figs. (JTJ) and (J2J)] and operators of Eq. © 
to tree level. Using dimensional regularization, where our subtraction scheme removes - + 
1 — 7e + log47r, we find the form factors have the form 

Gx{q 2 ) =Qx(l~ ) + J^jy 2 £ A i F H, q 2 ). (13) 

In this general expression, Qx is the meson charge, the sum on runs over all non-degenerate 
loop mesons of mass m^. The coefficients A* are products of charge and Clebsch-Gordan 
factors for the coupling of the loop meson to the external state meson X. These coefficients 
are listed in Table HI If a particular loop meson is not listed, then the values of A* are 
identically zero for all states X. Notice that in the isospin limit the charges qj and qi always 
enter in the combination qj + qi = —q r and for this reason do not explicitly appear in 
coefficients listed in the table. Lastly the non-analytic function F(m 2 ,q 2 ) above is defined 
to be 

1 [ 2 i m2 j_ a 2rf-<l 2 ' 
q log — + 4m y 7 



(m 2 ,q 2 ) = - 

D 



/i 2 \Am? 
where the auxiliary function T(a) is given by 



(14) 



I \l l — - + IS — 1 r 

T{a) = (a - 1) x j 1 - - log V a + -a - 2. (15) 

'1 - i + is + 1 6 



In the limit rrij m u , m r — > m s and using the physical quark charges, we recover the chiral 
perturbation theory result 

Using the values from the Table, the form factor of the 7r + , e.g., appears as 

G w+ (q 2 ) = jl - y^q 2 + j^jy [2F(m 2 u , q 2 ) + F(m 2 ru , q 2 )] J , (16) 
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TABLE I: Coefficients A% of loop mesons contributing to electromagnetic form factors. We list the 
coefficients only for X = ir + , K + , and K° states because the rest follow from charge conjugation: 
Af = -Af , A% + = -A$~, and A$° = -A~f for all <f>. And of course Af = A\ = 0. 



X 


ju ru 


i> 

js 


rs 


7T+ 


2(<? u - qd) q u ~ Qd 








K+ 


2q u + q r q u - q r 


-2q s - q r 


q r - q s 


K° 


2qd + q r qd- q r 


-2q s - q r 


q r - q s 



where the charge of the pion is Q n + = q u ~ Qd- The pion form factor appears rather special 
at one loop due to its independence from the charges of the sea quarks. 2 In fact, this 
independence from sea quark charges holds non-perturbatively and was shown using the 
behavior of lattice correlators under charge conjugation in ^3|. We demonstrate this more 
simply as follows. 

The electromagnetic current J M can be decomposed into isosinglet and isovector combi- 
nations, J M = + J 3 . Matrix elements of the former involve the charges of the sea quarks; 
while in the isospin limit, matrix elements of the latter are independent of these sea charges. 
Now consider the electromagnetic current matrix element of the neutral pion, 3 

(7r°(p')l J, \Ap)) = (Am 4 \Ap)) + (Am 4 k°(p)>. (17) 

This vanishes by charge conjugation invariance. Furthermore the isovector contribution 
vanishes, hence, so too must the isosinglet contribution. The matrix element of the isosinglet 
current is the same for all pions in the isospin limit. Thus as a consequence of isospin 
symmetry and charge conjugation invariance, we have 

<7T+(P')I ^ k + (P)> = <7T + (P')I Jl |7T+(P)>, (18) 

and hence the charged pion form factor is independent of the sea quark charges. 

As a result of this independence from sea quark charges, one can efficaciously ignore 
operator self-contractions and still determine the pion form factor [13(. Note that while 
the pion form factor is insensitive to contributions from closed quark loops with photon 
insertion, it is sensitive to sea quarks. As pointed out in [34], the analogous calculation 
in quenched chiral perturbation theory shows that there is no meson mass dependence at 
one-loop order. This result is also clear from the Table: only valence-sea loop mesons enter 
our one-loop expressions. Notice that the kaon form factors are not independent from the 



2 At one loop, we can see explicitly that the independence from sea quark charges arises from cancellations 
between isospin degenerate loop mesons. Away from the isospin limit the sea charges remain, e.g., in 
non-degenerate S'J7(4|2) the non-vanishing loop coefficients for the tt + are: = q u — qj, Af u = q u — qi, 

A jd = -(id - Qj), and = -(q d - qi). 

3 As with the proof in 13], our argument holds for charged SU(2) mesons of arbitrary spin and parity in 

the isospin limit, e.g. one can use the charge conjugation invariance of the p° to establish that the form 
factor of the p + is independent of the charges of the sea. 
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sea quark charges. This dependence only disappears in the SU(3) limit, which is badly 
violated in nature. 4 

Finally we derive expressions for the charge radii at one-loop order. In infinite volume, 
we can take q 2 <C m 2 to find that the meson charge radii are given by 

<4>=e*^-^£^(k>«5 + 1 )- (19) 



IV. LATTICE SPACING DEPENDENCE 

In this Section, we detail the modifications to our results at finite lattice spacing. We 
consider a mixed lattice action consisting of Ginsparg- Wilson valence quarks and staggered 
sea quarks. To address the effects of the lattice spacing, one formulates the continuum 
effective theory of the lattice action and then matches this effective theory onto a chiral 
perturbation theory. In this work we assume the natural hierarchy of scales 

1 

m q -C Aqcd < -, 
a 

and choose the power counting 

p 2 ~ ml ~ a 2 A 4 QCD . (20) 

The form factors can now be systematically calculated in the dual expansion in quark mass 
and lattice spacing. Such modifications to the electromagnetic form factors of mesons are 
rather simple. It was demonstrated in Ref. j4^| that no local a-dependent operators con- 
tribute to the form factors to 0(p 2 ) due to charge conservation. In our power counting, 
local corrections to the current in the a 2 chiral Lagrangian will contribute at 0(p 4 ) and will 
be competitive with two-loop effects. These local terms can thus be neglected here. The 
lattice spacing corrections to meson form factors then enter only through the dependence of 
the meson masses on the lattice spacing. 

With a mixed lattice action, there is no symmetry that relates the valence and sea sectors 
of the theor y. T he Symanzik Lagrangian thus contains dimension-six mixed field bilinears 

of the form jH, Q _ 

5C^ ~ {QTP V Q)(QTP S Q), 

where Py is a diagonal matrix that has unit entries corresponding to the valence sector, and 
Pg is a diagonal matrix that has unit entries corresponding to the sea sector (we have also 
implicitly enlarged the vector Q to transform under the continuum S77(15|3) that includes 
four tastes for each flavor of sea quark). Above V = 7^, 7^75 are the only Dirac matrices 
allowed by the chiral symmetry of the valence sector and the axial symmetry of the sea 
sector. The mapping of such mixed bilinears onto operators in the chiral perturbation 
theory at finite lattice spacing produces one new operator of the form [33| 
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5C = -ia 2 C mix str(T 3 £T 3 £t), (21) 

Indeed the argument presented in only applies for mesons consisting of degenerate flavors. In the 
non-degenerate case, the argument reduces to a demonstration that the sum of operator self-contractions 
in the meson and its charge conjugate vanish. 
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where T3 = Ps — Py- This operator contributes to masses of mesons formed from one valence 
and one sea quark. Such valence-sea mesons are not protected from additive mass renormal- 
ization due to the mixed action's symmetry. Thus to 0(p 2 ), the loop meson consisting of a 
sta gge red quark Q i} i.e. of flavor Q and taste i, and a valence quark Q' have masses given 

A a 2 
m QrQ' = -p( m Q + m Q') + -p C mix- (22) 

These masses are independent of the staggered quark taste and hence there is a four-fold taste 
degeneracy in the valence-sea meson loops. This four-fold degeneracy is exactly canceled by 
the j factor that must be inserted by hand to implement the fourth-root trick. Hence in a 
mixed action simulation, the meson form factors at 0(p 2 ) are given by Eq. ()13|) with the 
loop meson masses given in Eq. (|22|). 



V. VOLUME DEPENDENCE 

On the lattice, the available momentum modes are quantized and observables calculated 
thus inherit a dependence on the lattice volume. This dependence, which is inherently a long- 
distance effect, can be ascertained using chiral perturbation theory. We customarily choose 
a hypercubic box of three equal spatial dimensions L, and time dimension T, with T ^> L. 
With periodic boundary conditions on the quark fields, the available meson momenta have 
the form = (fc ,fc), with k = ^fn, and n represents a triplet of integers. Due to the 
assumed length of the time direction, we treat k as continuous. 

As spontaneous symmetry breaking does not occur in finite volumes, we must be careful 
also to specify that m n L > 1 so that mesonic zero modes do not become strongly coupled . 



This restriction ensures that the zero modes do not conspire to restore chiral symmetry. 
Provided this is the case, the standard p-counting of chiral perturbation theory remains 
intact and integrals over loop momenta can merely be replaced by corresponding sums over 
quantized momenta. This replacement leads to a dependence of physical quantities on the 
lattice size L. 

Calculation of the finite volume correction to our one-loop results for meson form factors 
is straightforward. Considering pionic matrix elements of the charge density 5 in a box, we 
have 

(Mp')\-MMp)) _ pn A _4a9 , 



Pa + Pi \ f 2 

eA x n 



1f 2 V> 



dx 



k + xq) 2 + m 2 , + x(l — x)q 2 Jk 2 + m 



,2 



The finite volume shift can then be obtained using the Poisson re-summation formula and 
the resulting sums can be cast into exponentially convergent forms, see, e.g. [45|. We find 



5 We do this only for simplicity. The form factors can be extracted from any component of the current, and 
this is the one commonly chosen in lattice QCD simulations. The volume effects from extracting form 
factors from the spatial components can similarly be calculated. 
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FIG. 3: Momentum transfer dependence of the pion form factor. The ratio of the pion form factor 
to its slope at zero momentum transfer is plotted versus q 2 for a range of pion masses. 

the finite volume shift to the meson form factors 



S L G x (q 2 ) = 




where we have defined 



l(q,A) = 




with $3(q, z) as the Jacobi elliptic theta function of the third kind. 

VI. EXTRAPOLATIONS OF THE CHARGE RADII 

We now consider each of the various extrapolations necessary to make physical predictions 
from lattice data for meson form factors. For simplicity we consider the pion form factor 
and this will enable us to compare with the lattice data in [lij . To perform our analysis, 
we fix the low-energy constant a 9 (/i = lGeV) = 0.0069 by using the one-loop result Eq. (JTHJ) 
along with the experimentally determined pion charge radius jlfjj]. 

A. Momentum Extrapolation 

First we shall deal with the limitation of discrete lattice momenta. Above we have 
determined the effect of the finite volume on the radii by using the allowed lattice momenta in 
a periodic box. Another consequence of periodic boundary conditions is that the momentum 
transfer is quantized; hence, the limit leading to Eq. (|19|) cannot be taken. Provided the 
lowest momentum transfers are in the chiral regime, the effective field theory can be used 
to perform a momentum extrapolation, cf. Eq. (|16j). It is questionable whether at current 
lattice volumes the minimum spatial momentum |q m in| ~ 0.5 GeV meets this restriction. For 
a relativistic object like the pion, a slight reduction is seen because q 2 < q 2 (unlike the 
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FIG. 4: Meson mass dependence of the pion charge radius. 

nucleon where q 2 ~ q 2 ). Furthermore, corrections to Eq. ()16|) from higher terms in the 
chiral expansion are of order g 2 /A 2 , for which the smallest available q 2 in Ref. yields 
~ 20% corrections. Thus we shall assume that the momentum transfer dependence of the 
lattice data at the minimal value of q 2 is captured by the effective field theory. 

With this assumption, let us investigate the momentum transfer dependence of the pion 
form factor at one-loop order. If the minimal momentum transfer q 2 is indeed small compared 
to m 2 , then the form factor will have linear behavior in q 2 and no momentum extrapolation 
is necessary to determine the slope near q 2 = 0. On the other hand, imagine that m 2 is 
small compared to q 2 . In this limit, the chiral logarithm dominates the form factor, which is 
multiplied by q 2 , and again the behavior is linear. In the intermediate region q 2 ~ 4m 2 the 
behavior of the function T\— g 2 /(4m 2 )] in Eq. (JT5j) becomes important. This is the region 
of parameters relevant for current lattice simulations. In Figure El we plot the function 
AG(q 2 ), defined by 

= ^gg^, (25) 

as a function of q 2 in order to see the deviation from linearity. The largest deviation from 
linearity over this range of momentum transfer is for the physical pion mass. For the larger 
pion masses employed on the lattice, the plot shows that for all practical purposes we can 
treat the form factor as linear in q 2 . For example, the data at the lowest g^ in and at lightest 
pion mass in Ref. 3|> the difference from linearity is less than one percent. Comparatively, 



the neglected higher-order terms in the chiral expansion are more than an order of magnitude 
larger (these terms, however, modify the momentum transfer dependence). 

B. Chiral Extrapolation 

Most striking is the well-known behavior of the charge radii near the chiral limit. This can 
be seen from the chiral logarithm in Eq. (|19p. For arbitrarily light pions the virtual cloud 
of pionic excitations in a meson extends arbitrarily far, and hence there is considerable 
variation of the charge radius as a function of the pion mass. One must be careful, however, 
not to assume that such variation with mass will necessarily be seen in lattice data. Looking 
at Eqs. (|19|) . and (|22j). we see that the valence-sea meson mass is what dominates the chiral 
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L [fm] 



FIG. 5: Volume dependence of the pion charge radius. The relative difference A of the charge 
radius in finite to infinite volume is plotted as a function of the lattice size L. The curve for 
m n = 0.14GeV terminates for m n L = 2 where pionic zero modes become important. 



behavior at one- loop order. Thus while the valence- valence and sea-sea meson masses are 
light, the additive mass renormalization allowed by the mixed- action symmetry breaking 
may actually push one away from the chiral regime. In Figure EJ we plot the meson mass 
dependence of the pion charge radius. Keep in mind that the meson mass which is relevant is 
the valence-sea mass. Ignoring other systematic errors, lattice calculations with valence-sea 
meson masses ~ 350 MeV will undershoot the pion charge radius by ~ 10%. 



C. Volume Extrapolation 

The virtual pion cloud will be affected by the boundary conditions imposed in the lattice 
simulation. Due to the chiral singularity in the one-loop contribution, one might also expect 
that the volume effects are substantial for light pions because these corrections stem from 
the long distance physics. Using Eq. (J23|) we can plot the volume effect as a function of L 
for various values of the (valence-sea) pion mass. In Figure 03 this is done for the quantity 
ArJ which is defined by 

where < > is the infinite volume radius, which is given in Eq. (jlDjl . and < >l= 
— ^^^lGtt(q 2 ) | 2=0 is the finite volume modification. We see that at the physical value 
of the pion mass, volume effects are substantial ~ 12% in a 3 f m box consistent with our 
intuition. However, the volume effects drop considerably for larger pion masses. For an 
m = 0.35 GeV valence-sea pion in a 2.5 fm box [which roughly corresponds to the values used 
in ignoring the effects of the mixed action in Eq. (J22)) ]. the effects of periodic boundary 
conditions in a finite box lead to a bigger positively charged pion, but by only a negligible 
0.5%. 
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D. Continuum Extrapolation 

Lastly we investigate the systematic error in the pion form factor data due to the finite 
size of the lattice spacing. The lattice spacing a enters our expressions only through the 
valence-sea meson masses, Eq. ()22|) . Ideally we would have data at multiple values of the 
lattice spacing and quark mass to enable a proper continuum and chiral extrapolation of the 
data. In Ref. [3], there is only one lattice spacing employed and only one valence- valence 
meson mass light enough to warrant a chiral analysis. As we have commented above, volume 
and momentum extrapolations produce negligible corrections compared to neglected higher- 
order chiral contributions. Thus we need only use the infinite volume expression in Eq. (jl6j) 
to make contact with the data. Fixing ag and / to their physical values, we can estimate 
the mixed action low-energy constant C m i x . 

Using the form factor data at the lowest value of g min = 0.18GeV 2 with the valence pion 
mass m n = 0.32 GeV and the sea pion mass rrijj = 0.35 GeV, we estimate C m i x = 0.0064 GeV 6 . 
We will not cite errors on this value because it is a rough estimate — only one datum is used 
to determine C m i x . Moreover as it enters our expressions logarithmically, any error in the 
form factor exponentiates into our estimate of C m i x . Using the statistical error bars from 
the lattice data for a high-low estimate, we find a rather wide range statistically allowed 
for C m i x from —0.0040 to 0.080 in units of GeV 6 . The range produced from the systematic 
uncertainty is comparable. 

VII. SUMMARY 

Above we have investigated various extrapolations necessary to connect lattice QCD 
data for meson form factors at small momentum transfer to the physical meson charge 
radii. We find that while volume effects are sizable at the physical pion mass, they are 
negligible for ~ 250 - 350 MeV in current lattice volumes. Provided the minimum lattice 
momentum g min < 300 MeV, there is no need for a momentum extrapolation: the small q 2 - 
dependence predicted from chiral perturbation theory is very linear. The corrections to this 
linearity arise from next-to-next-to-leading order contributions in the chiral expansion. At 
these pion masses and momenta, we thus conclude that systematic error is dominated by 
higher-order terms in the chiral expansion: ~ 10% from corrections that scale generically 
as m^/A 2 , and ~ 20% from terms of order <? 2 /A 2 with q 2 = g min of [l8j]. For the mixed 
action data in 0|, there are additional errors from lattice spacing artifacts. Demanding 
that their data be described within the effective field theory, we are able to estimate the low- 
energy constant C m i x that parametrizes the explicit breaking of the mixed action symmetry 
[577(613) -> 5C/(3|3) ® SU(3)} at 0{a 2 ). We obtained a rough value C mix « 0.006 GeV 6 , 
where there is sizable error in this value arising from the statistical error in the lattice data 
and systematic error due to neglecting of two-loop results (for both the quark mass and 
momentum dependence of the pion form factor). Because C m i x enters the expression for 
the form factor through a logarithm, a small ~ 10% statistical or systematic error in the 
lattice data results in a very large uncertainty for C m j x . Data at multiple lattice spacings 
or with varying quark mass would be ideal for a better determination of this low-energy 
constant. Nonetheless, the pion charge radius is a quantity that is sensitive to the fermion 
discretization, and can be used to understand the errors associated with the continuum 
extrapolation of mixed-action lattice data. 

The computation of meson form factors on the lattice is challenging. For the pion, 
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however, the absence of operator self-contractions in the isospin limit |13(, puts the charge 
radius in reach of current lattice technology. Further data at additional values of the valence 
quark mass and lattice spacing will thus allow one to predict the pion charge radius at about 
the ~ 10-20% level or better, without any model assumptions. 
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